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ABSTRACT 

We consider the impact of key astrophysical and measurement systematics on con- 
straints on dark energy and modifications to gravity on cosmic scales. We focus on 
upcoming photometric "Stage III" and "Stage IV" large scale structure surveys such 
as DES, SuMIRc, Euclid, LSST and WFIRST. We illustrate the different redshift 
dependencies of gravity modifications compared to intrinsic alignments, the main as- 
trophysical systematic. The way in which systematic uncertainties, such as galaxy bias 
and intrinsic alignments, are modelled can change dark energy equation of state and 
modified gravity figures of merit by a factor of four. The inclusion of cross-correlations 
of cosmic shear and galaxy position measurements helps reduce the loss of constrain- 
ing power from the lensing shear surveys. When forecasts for Planck CMB and Stage 
IV surveys are combined, constraints on the dark energy equation of state and modi- 
fied gravity model are recovered, relative to those from shear data with no systematic 
uncertainties, if fewer than 36 free parameters in total arc used to describe the galaxy 
bias and intrinsic alignment models as a function of scale and redshift. To facilitate 
future investigations, we also provide a fitting function for the matter power spectrum 
arising from the phenomenological modified gravity model we consider. 

Key words: cosmology: gravitational lensing: weak - dark energy - equation of state 
- cosmological parameters - large-scale structure of Universe 



1 INTRODUCTION 

Upcoming large scale structure surveys promise to provide 
deep and wide angular scale measurements of both the dis- 
tribution of luminous matter, through galaxy position sur- 
veys, and the total matter distribution, through weak lensing 
shear measurement. These include the ground-based Kilo- 
Degree Survey (KIDS), Pan-STARRS Q, Subaru Measure- 
ment of Images and Redshifts (SuMIRe) the Dark En- 
ergy Survey (DES) and the Large Synoptic Survey Tele- 
scope (LSST) 0, and prospective space-based Euclic0 and 
WFIRSlH surveys. The expectation is that these will sup- 
port a significant improvement in constraints on the effect 
of dark energy on the growth of cosmic structure. 

There has been an active discussion of the implications 
of such measurements for constraining cosmological models 
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that might predict growth histories that are distinct from 
those of the 'vanilla' ACDM cosmological model, in which 
the universe is populated by normal matter, cold dark mat- 
ter (CDM) and a cosmological constant (A) and an evo- 
lution determined by General Relativity (GR). Particular 
attention has been given to models arising from large scale 
modifications of gravity, such as those proposed as alterna- 
tives to expl ain the acce l eratin g expansion of the universe , 
for example iDvali et all d2000h; ISawicki fc Carroll! 
Zhand f2006ft;lAmarzguioui et all lpOOG); Bea n et ah] 
Movahed et all boqghlFang et al.l 
~Tcarvalho et al.l (120081 ) : ISchmidt et al 
lLombriser et al.l (I2009T) ; ICopeland et al 





iJain fc Khourv 
lHarko et al.l (|201ll ). 



(|2010h : ICapozziello et alj i|201lf ) 



Rather than considering a specific theory, one can 
consider a phenomenological parameterisation for how 
the gravitational metric perturbations are related to 
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Zhang et al. (2007); Linderl (20081); Jain & Zhane 
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(feOOSl); iBertschinger fc Zukinl J2008I); IZhao etaH (l2009bh; 
ISong fc Kovamal d2009h;lThomas et"aU (l2009h;ISong fc Pore] 



2009) 



il2009h; IZhao et al.1 (l2009bh : ISilvestri fc Trodden! 
Caldwell fc Kamionkowskil J2009I') ; Stril et al 
Zhao et al.l j2009al'); iGuzik et all TOld): iReves et al 



(2010) 



20101); iPogosian et al I j201Cf>: I Jennings et al J (l201ll) 
iMartinelli et al.l j201lF; iHwang et al.1 (l2010h; ISong et al 



IZhao 



2011 ) 



2011 ) 



d2010l);|Paniel et al. 1 d2010f> ; iBean fc Tangmatithaml d201oF 



et al.l d2010f): IPaniel fc Lindei 



Camera et all d201lf): Jain 



Clifton et al 



Linderl (|201lh ; iTsuiikawal (|2010l ); iBevnon et ~ 



20091 ). One parameterisation in particular is well suited to 
contrasting information about galaxy positions with weak 
lensing measurements. It specifies how the evolution of the 
two scalar metric perturbations in the conformal Newtonian 
gauge, the "Newtonian potentials", is distinct from that 
predicted by GR. While galaxy distributions are sensitive 
to just one of the potentials, lensing is sensitive to the sum 
of both. 

Cross-correlations between position and lensing observ- 
ables have been propo sed as a source of additional cosmo- 
logical inf o rmati on by iHu fc Jainl (|2004 ) and were used in 
iBernsteinl i2009h . It has been found that these observables 
could greatly incre ase constraining power for modified grav- 
ity (MG) theories (|Zhang et al.1 l2007t I2OO8I ) . however such 
cross-correlations also play an important role in mitigating 
critical sources of systematic error, in particular the effect 
of intrinsic align ments (IA) on the observed shear signal 
and galaxy bias (|joachimi fc Bridle! [20I0I : iKirk et al.ll2010l ; 
IZhandl201Ct ). 

However, since both an alternative theory of gravity 
and intrinsic alignments qualitatively alter the measured 
shear signal and its correlations with other observables, one 
might expect this to translate into a degeneracy between 
MG and IA parameters, a nd a deterior a tion o f cosmological 
constraints. For example, IZhang et al.l (|2007t ) highlight the 
use of galaxy-shear correlations (galaxy-galaxy lensing) to 
constrain modified gravi ty mod e ls in the absence of intrin- 
sic alignments, whereas IZhand (|2010T l proposes the use of 
the same correlations as a method for measuring intrinsic 
alignments in GR models. 

In this paper, we characterise and quantify the de- 
generacies between dark energy parameter measurements 
and weak lensing systematics, and assess how they af- 
fect our ability to obtain bounds on the nature of cos- 
mic ac celeration and gravity on cosmic scales. In a rel ated 
paper (| D. Kirk. I. Laszlo. S. Bridle and R. Beanll201ll ), wo 
discuss how surveys can be optimised in light of intrinsic 
alignment and galaxy bias systematics. 

Intrinsic alignments of galaxies arise from the fact that, 
as large scale structure is assembled, galaxies acquire an in- 
trinsic correlated ellipticity. Cosmic shear is the apparent co- 
herent alignment of galaxies due to gravitational distortions 
by intervening matter. Therefore, if intrinsic alignments are 
not accounted for then they can be falsely attributed to 
weak lensing and thus contaminate observations. Such in- 
trinsic correlations arise in two forms: an effect from physi- 
cally close galaxies forming or evolving in the same potential 
(termed II) and a correlation between the apparent shape 
of a distant galaxy lensed by a foreground mass and the 
intrinsic alignment of a foregr ound galaxy due to th e same 
foreground mass (termed GI) jHirata fc Seliakll2004l ). 



Including IA 
is well motivated 
tions reveals a need 
200d; 



effects in a cosmological analysis 
comparing simulations to observa- 



for these ef fects (ICroft fc Metzlerl 



Crittenden et al.l I2OO1I; iHeayens et al.1 |2000|; 



Hevmans et al.feood; iMandelbaum et al.1 20061; Hirata et al.l 



20071; lOkumura et all 120091; lOkumura fc Jind 12009^ 



Faltenbacher et al. 20091 ; Pereira fc Bryan 20ld ) and 



scale 


structure survevs (iBrown et al.l 


20021; iHirata et al. 


2004 


; Pereira & Kuhnl 20041; 


Mandelbaum et al. 


2005 


; Agustsson & Brainerd 20061; 


Hirata et al. 2007: 


Mandelbaum et al.l 20091; iBrainerd et all 20091; Siverd et al. 


2009; (Faltenbacher et al.l 120091; iLed 


2011; Joachimi et al. 


20ld; ted 120111; iBlazek et al.l l201ll; 


Hao et al.1 l201ll). At 



typical redshifts, the intrinsic alignment signal can be about 
ten per cent of the cosmic shear signal we want to measure. 
However, to measure the dark energy equation of state to 
percent level accuracy, we need to measure the cosmic shear 
signal itself to a similar accuracy, therefore incorrectly 
ignoring the intrinsic alignment signal can lead to biases of 
order tens of percent in the dark energy equation of state 
jBridle fc Kingj|2007f ). The difficulty in lensing then, even 
for a perfect instrument, is to separate out the true lensing 
signal (termed GG) from these intrinsic alignment terms. 

The II term can be removed by avoiding inclusion 
of physicall y close galaxy pairs as in iKing fc Schneider 
( | 2002 l , 120031 ) ; iHevmans fc Heavend (|2003l >; iTakada fc White 
l|2004r ). However, no similar simple treatment can be used 



for GI, since GI depends on line of sight interactions be- 
tween background alignments and forground lensing, and it 
is precisely the integrated lensing distortion along the line of 
sight that is relevant to the pure weak lensing signal. This 
leaves a contamination that must be handled in order for 
weak lensing to improve in quality even as in strumental er- 
rors are being beaten down (|Hirata et al.ll2007h . Recent work 
has shown that one can remove some IA contamination by 
"nulling" (|joachimi fc Schneider! l2010al lbl). This technique 
involves a transform of the shear data which downweights 
the impact of both II and GI terms. Nulling essentially re- 
moves information from the survey in order to produce un- 
biased, but necessarily less tight, constraints. An alternative 
approach is to model the IA and include additional param- 
eters to characterise the uncertainty in the model. 

There has been some effort centered on producing a 
robust theoretical model for IA formation and evolution. 



dCatelan et al. 2001 


; Hui & Zhang! |2002|; IHirata & Seliak! 


12004 Bridle & King 


120071; iLee et al.1 2008t). These models 



are based on the assumption that elliptical galaxies take 
shapes aligned with those of their parent dark matter haloes, 
which themselves respond linearly to the background gravi- 
tational tidal field. Spiral galaxies have shapes produced by 
torques acting on the angular momentum axis of the galax- 
ies and are shown not to produce IAs to first order. The LA 
model is a simple model which generates results consistent 
with GI measurements (|joachimi et al.l l201ll ; IBlazek et al.l 

Although the linear alignment model serves as a suit- 
able starting point for generating agreement with observa- 
tions, there are several complicating factors. IAs are likely to 
be modified through merger and accretion events as struc- 
ture forms, and also impacted by additional factors such as 



© 2009 RAS, MNRAS OOO.riHTsl 



Disentangling dark energy and cosmic tests of gravity from weak lensing systematics 3 



color, redshift, galax y type, and luminosity dependence and 
one-halo corrections dHirata et aL 2007 Faltenbacher et al.l 



Schn eider 
201 il l 



20091; iMandelbaum et al.ll2009l; lHao et al.1 l201ll; lLedl2oilI 



Bridle! 1201(3 : iKirk et all l20ld ; Ijoachimi et all 



Rather than attempt to model each of these complex 
effects directly, we allow for added freedom in the LA model 
by introducing a gridded bias scheme, gridded in scale and 
redshift space, and marginalise over the additional parame- 
ters. Including such a scheme to account for uncertainties in 
the IA amplitude, variation with galaxy type and evolution 
effects, wi l l imp a ct constraints on cosmological parameters 
iBernsteinl d2009h; iBridle fc Kind (|2007f l; iKitching fc Taylorl 



(20 1 lis 



Joachimi et all (|201lh . 



A primary goal of future surveys is to beat down instru- 
mental systematic errors to improve cosmological constrain- 
ing power. In extracting cosmological constraints, however, 
one has to also disentangle astrophysical systematic effects. 
Along with the effect of baryons on small scales, our ability 
to model the growth of non-linearities in large scale struc- 
ture and the contamination due to intrinsic alignments are, 
perhaps, the main astrophysical systematic challenges facing 
weak lensing surveys. The formation, evolution as a function 
of redshift and scale, and variation by galaxy type of IAs are 
all relativity unknown. Realistically one might expect that 
one will have to extract an understanding of dark energy 
and IAs simultaneously from lensing and galaxy position 
correlations, this is particularly so in the context of modi- 
fied gravity, which affects both the cosmic homogenous ex- 
pansion and inhomogeneous perturbation growth histories. 
Here we present a detailed analysis of this task at the level 
of a Fisher matrix analysis. 

The paper is structured as follows: in section [5] we out- 
line the formalism used to describe the cosmological model, 
including dark energy parameters to describe modifications 
to ACDM in how the homogeneous expansion and growth 
of inhomogeneities evolve. We describe how the large scale 
structure and CMB surveys are modeled and how we treat 
uncertainties in galaxy bias and intrinsic alignments. In sec- 
tion [3] we present the results of the Fisher matrix investiga- 
tion and summarise our findings in [4] In an Appendix we 
introduce a fitting function that allows others to generate 
weak lensing and galaxy position correlations for the modi- 
fied gravity model we consider here. 



2 FORMALISM 

In this section we describe the details of our analytical ap- 
proach: the cosmological model we use to parameterise dark 
energy and modifications to gravity is outlined in sect ion [2~T1 
sections 12.21 and |2 . 31 respectively describe the statistical cor- 
relations for large scale structure and CMB observables and 
the survey specifications we assume. In 12.41 and 12.51 we de- 
scribe how systematic uncertainties in the galaxy bias model 
and intrinsic alignments are included in the analysis. 



2.1 Cosmological model 

Theories that suggest a cosmic scale modification to gravity 
can produce deviations from ACDM in both the homoge- 
neous expansion history and the growth of inhomogeneities. 



We model alterations from ACDM in the expansion history 
using an effective equation of state for the additional physics 
producing the cosmic acceleration 



w = wo + w a (1 — a), 



(1) 



where a is the expansion factor, with a = 1 today, and wo 
and w a describe the effective equation of state today and 
its derivative with respect to a respectively (wo — — 1 and 
w a = recover a cosmological constant, A.). 

The Friedmann equation relates the Hubble expansion 
rate, H(a), to the cosmic energy density, p(a) 



H 2 (a) = 



8tvG 



Hi 



— =- H r + n A a 

a- 1 a 4 



3(l+WQ-\-w a ) ^3w a (1 — a) 



(2) 



(3) 



where pi is the homogeneous (background) density of com- 
ponent i, Ho is Hubble's constant, Q m and fi 7 are the frac- 
tional energy densities today in non-relativistic and rela- 
tivistic matter respectively. 

To describe the growth of inhomogeneities we use 
the conformal Newtonia n gauge using the notation of 
iMa fc Bertschingerl (|l995l ) 



-a (r)[l + 2i>(x,r)]dT + a (r)[l - 2<p(x, T)]dx 2 (4) 



where r is conformal time, Xi are comoving coordinates and 
ip and <j> are the Newtonian potentials. 

We consider a phenomenological parameterisation for 
cosmic scale deviations from GR that employs two functions 
to modify the perturbed Einstein equations. Such modifica- 
tions have been widely discussed in the literature, both be- 
cause they are phenomenologically simple, and also because 
they can be mapped onto predictions for modifications de- 
rived from scalar tensor and higher dimensional theories of 
gravity. 

Unlike with the equation of state, notation for these 
modifications varies widely. Here we write the modified Ein- 
stein equations as 



k 2 4> = -4-irGQa 2 P»A t 

i 

tP~R4> = -127rGQa 2 ^pi(l+w)^. 



(5) 
(6) 



where the first, the Poisson equation, comes from a combina- 
tion of the time-time and time-space Einstein equations, and 
the second is the anisotropic space-space equation. Here k is 
the comoving wavenumber, Aj is the perturbation of compo- 
nent i in the component's rest-frame (the frame in which it 
has zero peculiar velocity) and Oi is its anisotropic stress. At 
is a gauge invariant perturbation, related to a perturbation 
in a general frame, 5t, in which the component's peculiar ve- 
locity is, Vi, by Ai = Si + 3H(1 + Wi)vi/k. In the discussion 
below, we focus on density perturbations in CDM, denoted 
by A c . 

In General Relativity the two functions Q — R = 1. The 
function Q(k,a) describes the relation between space-space 
gravitational potential perturbation and a matter overden- 
sity through the Poisson equation, creating an effective 
gravitational constant G e //(fc,a) = Q(k,a)G. The second 
function R(k, a) modifies the relationship between the two 
Newtonian potentials. In GR, an inequality between these 
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can only be generated by the presence of relativistic mat- 
ter, through the presence of anisotropic shear stress. Shear 
stresses are rapidly suppressed in non-relativistic matter, 
and would be extraordinarily difficult to sustain in a fluid 
with a negative equation of state. As such, R ^ 1 could be 
viewed as a potential smoking gun signal of a modification 
to GR being present. We assume that the fluid equations for 
normal matter are unaltered and assume that dark energy, 
as a modification of gravity rather than a fluid, does not 
cluster. 

The growth of overdensities on large and small 
scales have different depe ndencies on Q and R 
iBean fc Tangmatithaml (|2010h . On small scales, they 
are purely determined by the peculiar motion during infall, 
ip, proportional to QR, 



A c + HA C 



n m (a)QRA c w 0, 



(7) 



where "H = dlna/dr — aH . For late-time, large scale be- 
haviour the degeneracy between Q and R is not present and 
CDM density perturbation evolution is governed by 



A c + HA C - 3(0 + n{2(j) + </>) + (2H + H )4>) « 0. 



(8) 



We assume that time evolution of the effects of modifi- 
cations to GR on the growth of inhomogeneities would vary 
in concert with alterations to the background expansion, and 
consider modifications of the form, 



Q(a) = l + (Q -l)a s 
R{a) = l + (Ro- l)a s 



(9) 
(10) 



For our analysis we fix s — 3 to allow modifications to the 
growth of structure to evolve at a comparable rate to the 
onset of an accelerative component in the homogeneous ex- 
pansion history. We omit any scale dependence in our mod- 
ifications. Two conditions, Q(a) > and R(a) > — 1 are 
imposed to ensure that overdensities remain gravitationally 
attractive and that light is bent towards the lens. 



2.2 Observables 

We are interested in how well measurements of weak lensing 
shear distortions and galaxy counts, from future surveys, 
will be able to constrain deviations from ACDM. Combining 
these observables is key to detecting changes in the growth 
history of the universe as they depend differently on the 
Newtonian potentials. 

We characterise each observable by their 2D angular 
power spectrum, Ce, for auto- and cross-correlations be- 
tween each observable in a given redshift bin. For two fields 
X and Y, Ce, under the Limber approximation, is given by 



C 



° ^Wx(x)w Y (x)Sx(kt,x)SY(kt, x ) (ii) 

X 2 

where kt = i/x> an d X, Y={8, G} for mass and lensing shear 
fields respectively. Wx and Sx are the window function and 
source function associated with the field X, respectively. 

The source functions for the mass distribution and weak 
lensing shear are 



Si = A c , 
So = -yW> + V0- 



(12) 
(13) 



where we have assumed that the density perturbation for 
matter is equivalent to that for CDM 

The power spectrum for the correlation between fields 
X and Y is related to the source functions by, Pxy = 
(SxSy), allowing one to rewrite the angular correlation in 
a common form: 



C 



-4Wx(x)W Y (x)PxY(h, X )- 



(14) 



The galaxy and lensing window functions are dependent 
on the normalised distribution of galaxy number density in 
each redshift bin i for the relevant survey, rii(\). We assume 
the galaxies are distributed according to (|Smail et al.lll994h 
with 



n(z) oc z exp ( 

20 



3/2 



(15) 



We break up the galaxy distribution into N p h photometric 
redshift bins, divided so that they each contain an equal 
number of galaxies. The photometric redshifts are measured 
with accuracy cr(z) = o" z o(l + z) and could have a potential 
systematic offset, which we model as a constant within each 
redshift bin, Azi. The observed distribution of galaxies in 
bin i is given by 

rii(z) = —7^- [erf(«i — z + Azi) — erf(zj_i — z + Azj-i)] 
The galaxy and lensing window functions are then given by 



WUx) = n % (z) 

W b (x) ^ £^*W)^y 



x) 



(17) 
(18) 



where r(x) is the comoving angular diameter distance to co- 
moving distance x, an d hi is the normalised number density, 



hi(z) = 



m(z) 



(19) 



In relating the observed correlations of galaxies to the 
dark matter correlation functions above, we must account 
for a bias between dark and luminous matter. This bias is 
dependent on the galaxy type, redshift and environment in a 
way that is poorly understood. We allow for this uncertainty 
by introducing a redshift and scale dependent bias parame- 
ter, b g , to relate the auto-correlation of the galaxies to the 
autocorrelation of the mass and an independent, correlation 
parameter, r g , to describe the bias in the cross-correlation 
of luminous matter and the mass. The correlation param- 
eter, r g , would be equal to unity if the galaxies trace the 
mass, and less than unity if there i s some stochasticity in 
the galaxy formation process e.g. see lDekel fc Lahavl (1 19991 ) 
for a review. 

The galaxy position auto-correlation and galaxy 
position-shear cross-correlation are then related to the un- 
derlying mass and shear observables by 



p ss( fc .x) = b g (k,x)Pss{k,x), 
Pgc(k,x) = b g (k,x)r g {k,x)Psa(k,x)- 



(20) 
(21) 



where g and 5, denote galaxy and underlying mass observ- 
ables respectively. 

Note that in modified gravity models the lensing source 
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Survey Parameters 


Stage III 


Stage IV 


Arca(sq. deg.) 


5000 


20000 


V2z 


0.8 


0.9 


lit 1 It 


0.001 


0.001 


Zmax 


3 


3 


N g 


10 


35 


N ph 


5 


10 


&zO 


0.07 


0.05 


^frms 


0.23 


0.35 



Table 1. Summary of the photometric large scale structure sur- 
vey specifications assumed for the Stage III and Stage IV survey: 
survey area; median survey redshift, \/2zq; minimum and maxi- 
mum redshifts observed, z m i n and z m ax', number of galaxies, per 
square arcminute, N g ; number of photometric redshift bins, N p h\ 
standard photometric redshift measurement error at z = 0, cr z o, 
and the r.m.s. shear measurement error, ^rms • 



!/(GHz) 


100 


143 


217 


f sky 


0.8 


0.8 


0.8 


0FWHM(a?c min) 


10.7 


8.0 


5.5 


a T {fiK) 


5.4 


6.0 


13.1 


<7E (/iK) 




11.4 


26.7 



Table 2. CMB survey specifications for a Planck-like survey. We 
model this on the temperature, T, and E-mode polarisation spec- 
ifications from three lowest frequency bands for the Planck HFI 
instrument. 



term and the mass source term must be different 

~Q(x)(R(x) + i) 



Psc(k,x) = 
PGG(k, X ) = 



Q(x)(R(x) + 1) 



Pss(k,x), 



Pss(k, X ). 



(22) 
(23) 



The growth of the dimensionless power spectrum Pss is it- 
self dependent on modified gravity parameters Q and R, as 
summarised by (0 and JS}. 

To obtain the lensing and galaxy position correlations in 
the modified gravity scenarios we integrate the f ull equations 
of motion using a modified version of CAMB (E ewis et al.l 
l2000h . 

To support other researchers investigating the role of 
modified gravity models on large scale structure observa- 
tions, without having to integrate the full perturbation equa- 
tions, we provide a fitting function in the Appendix for 
the ratio, r/; t (fc, z), between a fiducial ACDM linear matter 
power spectrum, Pss,ACDAi(k, z) and the one for a modified 
gravity model described in l2.ll parameterised by Qo, Ro and 



rfit(k,z;Q ,Ro,s) 



Pss,fit{k, z; Qo, Ro, s) 



:,ACDNl(k, z) 



(24) 



2.3 Survey specifications 



We consider the impact of including IAs on cosmological 
constraints for a near- term Dark Energy Task Force (DETF) 
lAlbrecht et~al1 (|2006h Stage III survey, such as DES or 
SuMIRe, and a longer-term Stage IV survey, such as Eu- 
clid, LSST or WFIRST. 



The noise for each survey is modeled as statistical errors 
given by 



N Htj _ r _7rms 



2m 



N 



0. 



(25) 

(26) 
(27) 



where 7 rms is the root mean square uncertainty in the shear 
measurement of the galaxies and rij is number of galaxies 
per steradian in j th photometric redshift bin so = N g . 

The survey specifications assumed in our analysis for 
the Stage III and IV surveys are given in Table [T] 

We include complementary constraints from tempera- 
ture (T) and E-mode polarisation (E) measurements from 
a Planck-like CMB survey up to I = 3000. As summarised 
in Table [2] we model this by considering the three lowest 
frequency bands of the Planck HFI instrument, three chan- 
nels for temperature data and 2 for E mode polarisation, as 
described in the Planck Bluebook Q. We assume each fre- 
quency channel has Gaussian beams of width Qfwhm and 
error in X = T, E of ax , so that the noise in channel c is 
given by 



N, 



XX, c 



5c „„„„ \2/(f+l)eFWffM,c/8M2) 
'F\\ HM.cl £ , 



and over all channels, 



n; 



£0 



CT 1 



(28) 



(29) 



2.4 Intrinsic Alignments 

Cosmic shear describes the distortion of the image of a dis- 
tant galaxy due to the bending of light from that galaxy 
by gravity as it passes massive large-scale structure. For a 
galaxy in the i th photo-z bin, the observed ellipticity, e, of 
the galaxy can be written as a sum of three independent 
contributions: the cosmic shear 'ya, the intrinsic, non-lensed 
shape of the galaxy, 71, and apparent ellipticity introduced 
through instrumental and foreground noise, e rrl£ j, 



e\6) 



7g(0)+t1(0) 



(30) 



The cosmic shear signal 7c is very small, and we cannot 
measure directly the intrinsic shear of any individual galaxy. 
To recover the cosmic shear, therefore, one averages over a 
number of galaxies on a small patch on a sky. Assuming 
that their intrinsic ellipticities are distributed randomly, and 
that their light passes by similar large scale structure, the 
intrinsic ellipticities cancel in the two-point function, and 
we are left with the cosmic shear signal. 

In reality, the assumption that intrinsic ellipticities are 
randomly distributed on the sky is inaccurate. There are 
two strains of intrinsic alignment of galaxy ellipticities, both 
arising from the same physics of galaxy formation. 

The measured weak lensing signal reflects a correlation 
in shapes arising from distant galaxies passing near the same 
foreground gravitational lens. However, if the background 



7 www.rssd.esa.int/ S A/ PLANCK /docs /Bluebook 
SCI (2005) l.V2.pdf 
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images are already correlated, this boosts the measured sig- 
nal and gives rise to a systematic deviation in the observed 
shear. 

In section 12.4.1 we describe the linear a lignment 
model jCatelan et all 1200 ll ; Iffirata fc Seliakl |2004 ) we use 
as our basis to describe intrinsic alignment contribu- 
tions. In this model, galaxies quickly align along the cur- 
vature of the gravitational potential. While this model 
does provide a reasonable starting point for an overall 
agreement with observations, and thus a suitable base- 
line for use in our statistical analysis, it neglects known 
complicating factors of dependency on galaxy luminosity, 
galaxy type and redshift, and the effect of po st-processing 
of IA in mergers and accre t ion events (Irlirata et al.l 
20071; iFaltenbacher et "all 120091; iMandelbaum et al.l \20m 



Hao et all l201ll; ISchneider fc Bridle! |2010| ; Ijoachimi et al.1 
201ll ; lLedl201ll ; lBlazek et alj|201ll )."ln section [2A2l we there- 
fore extend beyond the LA model, in an attempt to allow for 
these effects, through the inclusion of additional non-linear 
corrections, and scale and redshift dependent bias terms. 



2.4-1 The Linear Alignment Model 

The Linear Alignm e nt (LA) mode l intro duced in 
ICatelan et"aH (120011 ); lllirata fc Seliakl |2004 ) assumes 
that galaxies would align with the stretching axis of the 
potential in which they form, so that the intrinsic shear is 
assumed to be 



(1,2) 

7} 



Cf 
4ttG 



(v 2 



vj,v t 



(31) 



ip(zf) is the smoothed gravitational potential field sourc- 
ing the shear alignments at a primordial redshift, z/,when 
galaxy formation occured and Cf is a normalisation de- 
termined by matching the model to observations. We take 
Zf =50 here, however we expect Zf to be well within the 
matter dominated era, when the potential ip remains roughly 
constant in time and gravity can be well-described by GR, so 
that the analysis should be largely insensitive to the precise 
value of Zf assumed. 

The presence of intrinsic alignments alters the observed 
shear correlations. We denote the measured galaxy shear and 
position observables by e and n, to distinguish them from 
the underlying theoretical IA-free shear and position vari- 
ables, G and g respectively. The observed angular correla- 
tion functions for galaxy shear and position, when intrinsic 
alignments are included, are given by 

CT = Cf G + Cf 1 + Cl 1 , (32) 
CT = Cf + Cf (33) 

ignoring magnification terms (|joachimi fc Bridiell2010l ). 

The Intrinsic-Intrinsic (II) alignment correlation ap- 
plies to physically close galaxies. These form in the same 
large-scale gravitational potential and their intrinsic ellip- 
ticities tend to align with the field lines of that potential. 
When these galaxies are observed on the sky they will tend 
to point in the same direction. This alignment produces a 
spurious correlation which adds to the observed cosmic shear 
signal. 

Since the II correlation is greatest for closely positioned 
galaxies, there is a related correlation between the position 



and the intrinsic alignment of a pair of physically close galax- 
ies, this is the galaxy position-intrinsic alignment (gl) cor- 
relation. 

Somewhat more subtle is the Gravitational-Intrinsic 
alignment (CI) correlation which applies to galaxies close 
on the sky but separated in redshift. The intrinsic ellip- 
ticity of the foreground galaxy will tend to align with the 
nearby gravitational potential which, in turn, is responsi- 
ble for gravitationally lensing the background galaxy. This 
tends to produce an anti-correlation which subtracts from 
the observed cosmic shear signal as we observe the galaxies 
oriented orthogonally. 

We can write angular correlation functions involving in- 
trinsic alignments using the formalism in (|11[) . denning a 
source term, Si, and window function, Wi, for the intrinsic 
alignments 



Si(k, X ) 

wi( X ) 



Cf 
4tyG 



k ip(k,z f ), 



(34) 
(35) 



The angular correlation is often written in terms of the linear 
power spectra, P X y = (SxSy), 



Cl*{l) 



%WaAx)Mx)PGi(k, X ), 
X 

%n l ( X )n J (x)Pii(k,x), 
X 2 

X 2 



(36) 



;n i (x)h j {x)P s i{k,x)- 



One can write the correlations in (|36|l in terms of the matter 
power spectrum at Zf, 



Pai(k,z) 
P u (k,z) 
P gI (k,z) 



Q(k,z)[l + R(k,z)] D(k,z) 
2 D(k, Zf \ 

C 2 f p 2 m (l + Zf) 2 Ps S (k,Zf), 

D(k z) 

D(k,z f ) 



Cfp m {l + Zf)P 6S (k, Zf), 

(37) 



where p m is the mean matter density today and D(k, z) is 
the linear growth factor for CDM perturbations, 



Ac(M) = ^4a c (M/) 
D(k,Zf) 



(38) 



In GR the linear growth factor is scale independent. If 
gravity is modified on cosmic scales, however, it can be scale 
dependent, and sensitive to the functions Q and R. 

How does one normalise t he IA correlations, and the 
constant Cf ? iBridle fc Kmel (|2007l) provide a numerical 
value for C/ 7 , by comparing with Hirata fc Seliakl (|2004T l 
who used observations of low redshift galaxies today. Rather 
than normalising the IA source function at Zf, they nor- 
malise it, with constant Ci, relative to the gravitational po- 
tential today assuming a ACDM cosmology, ?/>a, 



&(fc,z) = 
They find Ci 



-^k 2 Mk,o). 



(39) 



5 x 10- 14 (/i 2 M slm /Mpc- 3 )- 



8.25fe 2 x 10 4 MjPC 2 . We extrapolate the normalisation of 
ijBridle fc Kindl2007h at z = to z = Zf by assuming a 
ACDM growth factor, D A (z) (|Carroll et al.lll992l) . The early 
and late time normalisations, Cf from C\, are then related 
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by 
Cf 



-Da(O) 



D A (z f )(l + z f Y 



(40) 



2.4-2 Generalising the IA model 

Though the correlations in galaxy orientation might be 
formed in the era of galaxy formation, the intrinsic align- 
ments we observe will invariably be sensitive to the evolution 
of those galaxies, and the galaxy environment. 

One factor that is not included in the LA model is the 
impact of non-linear clustering of galaxies on the distribu- 
tion of the galaxies sourcing the intrinsic alignments. A non- 
linear alignment (NLA) model was introduced as an ad-hoc 
way to incorporate non-linear clustering into the LA model 
(|Hirata et al.ll2007l ; iBridle fc Kindl2007f ). This replaces the 
linear power spectrum in the intrinsic alignment angular cor- 
relation expression (|37[) with the non-linear power spectrum 
based on the fittin g function derived from the halo model of 
ISmith et all (|2003h . 

While this is somewhat adhoc, it was found to give 
a qualitatively similar result to the more motivated halo 
model of intrinsic alignments in ISchneider fe Bridle! (|2010| ) 
and it has been sh own, in for example Hirata et al.l (12007 ' 



iBlazek et al.l (|201ll ). that it gives a more consistent fit to the 
data than the linear alignment model. 

We extend this approach to take into account the mod- 
ified growth history, if gravity deviates from GR. We model 
the effect of non-linear clustering by boosting the IA source 
function, Sj fLA relative to that predicted by linear align- 
ments, 



S? LA (k,x) = Sf A (k, X ) 



Pss(k,x) nonHn 



Pss(k,xY 



(41) 



here Pjj n is the linear matter power spectrum predicted 
modified gravity model, and Pss" m is the non-linear spec- 
trum after a correction is applied to the pow er spectrum 
using the Smith et. al. halo fitting function (|Smith et al.l 
l2003h . The use of the Smith et al. fitting function, for a mod- 
ified expansion history, is equivalent to assuming that non- 
linear collapse in the modified gravity theories follows the 
Zel'dovich approximation; this was shown to be reasonable 
if the phenenomolo gical modifications in ([5]) and Q hold 
to nonlinear scales (IStabenau fe JainllioO a: iLaszlo fc Bear] 
120081 ). We briefly discuss the motivation and possible impact 
of deviations from this assumption later in the analysis. 

In figure [1] we show how each intrinsic alignment con- 
tribution to the observed correlations varies as a function 
of redshift for a fix ed multipole, £ = 1 000. We refer the 
interested reader to Ijoachimi fc Bridlel l)20ld ) for a figure 
detailing the variation in lensing, galaxy and IA correla- 
tions across all multipoles and redshift bins. The same-bin 
correlations vary monotonically as a function of redshift (de- 
noted by the index of the photometric redshift bin), while 
the cross-bin correlations involving galaxy positions dramat- 
ically fall off as the photometric redshift bins become more 
separated. By contrast, the cross-bin correlations for GG 
and GI, can remain significant even in cross-correlations be- 
tween widely separated bins because of the broad redshift 
kernel for the lensing window function. 

Figure [2] compares the amplitude and variation in shear 



and galaxy correlations when, separately, I As are added and 
when the modified gravity theory is allowed. While the two 
effects can be of comparable amplitude, their distinct red- 
shift dependencies, if the intrinsic alignment model is per- 
fectly known, could assist in disentangling them. 

Intrinsic alignments are unlikely, however, to be de- 
scribed fully by the linear alignment model; they will de- 
pend on the details of galaxy formation within dark mat- 
ter halos, and baryonic physics within galaxies, with com- 
plexity beyond this model. For example, it is known that 
the IA sig n al depends strongly on galaxy type and color, 
(|Led 12011 ; IBlazek et alj l201ll ); spirals are supported by 
angular momentum, and so more likely subject to tidal 
torquing of the angular momentum vector, while ellip- 
tical galaxies are better described by linear alignments. 
This bifurcation translates into colour dependence, as spi- 
rals are blue while ellipticals are older and r edder, and is 
noted in surveys which sp l it samples by colour (iHirata et al.l 
120071 ; IFaltenbacher et al.l 120091 ; IMandelbaum et all I2009T T 
Redshift and luminosity dependences, and one-halo cor- 
rections at smaller scale s also exist, as in (IHirata et all 



20071; IFaltenbacher et~ai1 |2009|; IMandelbaum et al 



2009 



Hao et all l201ll; ISchneider fc Bridlel |2010|; Ijoachimi et all 



20111) . 



Rather than attempt to incorporate these numerous ef- 
fects by direct modelling, we choose to allow a scale and 
redshift dependent bias factor to parameterise our ignorance 
and marginalise over the bias parameter in a redshift and 
scale gridding. We introduce two additional bias parameters 
into the IA correlations, bi and 77, analogous to the galaxy 
bias b g and r g , to reflect our uncertainty in the bias model: 



Pn = 



Pgi = biriP™ 



(42) 
(43) 



where the fiducial functions are given by the LA or NLA 
mode. 



2.5 Modeling galaxy bias and IA amplitudes 

We consider two scenarios for galaxy and IA bias. In our 
simple model all bias parameters are k and z independent, 
i.e. b g , bi, r g and ri are each a single constant free parameter 
determining the amplit ude of the bias. In the mo re realistic 
scenario, motivated by Ijoachimi fc Bridlel (|2010h each bias 
coefficient, Bx £ {b g ,r g ,bi,ri}, is interpolated from a Nk X 
N z grid of values logarithmically spaced in k and z, B x , 
respectively, 

B x (k,a) = (1-A fc ) [{l-A z )B x +A z B x +1 } 



[< 



+A fe \ (1-A,)B 



+ A Z B 



i+lj + l 



for ki < k < ki+i and Zj < z < Zj+i, with 

Ak _ ln(fc/fc t ) 
\n(k i+1 /ki) 

ln[(l + *)/(l + Zj)] 
ln[(l + ^ +1 )/(l + 

and modulated by a free constant amplitude parameter. 
We choose k m in = 10~ 3 Afpc _1 and k ma x = 30Mpc _1 
for the gridding, and assume Bx(k < k m i n ) = 1 and 
B x (k > k m ax) = B 



(44) 

(45) 
(46) 



N k N z 



at all times. We consider scenarios 
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123456789 10 123456789 10 

N ; N, 

Figure 1. A comparison of the intrinsic alignment and cosmological contributions, assuming a fiducial ACDM cosmology, to the shear- 
shear [left panel] and position-shear correlations [right panel] as a function of photometric rcdshift bin, jVj, for the Stage IV survey 
specification for a single multipole, I = 1000. Same-bin Hi' [full lines] and cross-bin correlations with the 5 th , central, rcdshift bin '5i' 
[dashed] are shown for the cosmological correlations GG, gg and gG, and the intrinsic aligment correlations 77, GI, and gl. For the gG, 
gl and GI correlations we plot the larger of i5 and 5i correlations in each case. 




io~ a l , , , , , , , , ] I , , , , , , , , U 

123456789 10 123456789 10 

N; N; 



Figure 2. The difference in shear-shear [left panel] and position-shear correlations [right panel] in which IAs are included [blue lines] 
or using a modified gravity model (with no IAs) [red lines] in comparison to a fiducial model, in which no IAs are included and GR 
is assumed. The modified gravity model shown has Qq = 1.05 and Ro = 1. As in figure [T] we show correlations for i = 1000 in each 
photometric rcdshift bin, JVj, for the Stage IV specification. 



in which 1 ^ Nk = N z ^5. Each of the Nk x N z grid nodes 
is a freely varying parameter. This means that, in our more 
sophisticated model, there are AN^N Z nuisance parameters 
when all correlations and cross-correlations are included. 

A multi-bin marginalisation over bias parameters is ar- 
guably conservative, however we believe it reasonably re- 
flects the current uncertainties in the bias and IA models. 
It was inspired by the work of |Bernstei"r] (|2009l ) which was 
used in lAlbrecht et all (|2009h who bin the biases in redshift 
and multipole bins, rather than redshift and wavenumber as 
we do here. 

For observables involving the galaxy position correla- 
tions we truncate the maximum i us e d as a function of 
redshi ft bin, as per lRassat et al.l (|200et ); Ijoachimi Sz Bridle] 
(|20ld ) to remove poorly understood biasing on non-linear 
scales from the likelihood calculation. We introduce a max- 
imum wavenumber k max for a photometric redshift bin i, 
and neglect all ti > k max x(zi). 

Our analysis spans from the optimistic to conservative 



scenarios. Optimisitically one might assume one can extract, 
and therefore exclude, IAs with perfect precision and can 
model the galaxy bias as scale and redshift independent. A 
conservative perspective would be to represent our ignorance 
in IA and galaxy bias modeling with 100 marginalised pa- 
rameters (Nk = N z = 5). Where, within this range, the 
realistic range will finally fall will depend on progress in un- 
derstanding IAs potentially through the use of complemen- 
tary spectroscopic redshift surveys and the development of 
galaxy training sets, or preferential selection of galaxy sub- 
groups in which intrinsic alignments are less pronounced, 
and galaxy bias is well-understood. 



3 ANALYSIS 

For our analysis of the impact of systematics on dark en- 
ergy constraints, we consider constraints on 10 cosmological 
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parameters: 



Qo{l + Ro) 



P = {Qa, WO, W a , Qo, V " UV ^ — , fib/l 2 , fi m /l 2 , T re i on , 

n s ,ln(10 10 yl s )} (47) 

where T re ion is the optical depth to the epoch of reionisation, 
and n s and y4 s are the spectral index and normalisation of 
the primordial spectrum of curvature perturbations, with 
pivot scale k — 0. 05 M pc~ 1 . We choose fiducial values for 
these parameters assuming ACDM, and consis tent with a 
WMAP7 bestfit cosmology jLarson et al.ll2010l ). 

As can be seen in !2.2l a primary constraint on modified 
gravity parameters from weak lensing is the combination 
from Qo(l + Ro)/2, rather than Ro- We therefore use Qo 
and Qo(l + -Ro)/2 as variables in the Fisher analysis, and 
take GR with Q — R = 1 as the fiducial model. 

Unless stated otherwise, we consider a conservative sce- 
nario for astrophysical systematics, and marginalise over 
the galaxy and intrinsic alignment biases b g ,r g ,bi,ri with 
Nk — N~ = 5 bins. We assume b g — r g = br = Ti = 1 for 
the fiducial model. 

Our fiducial scenario involves a 5 x 5 grid of 4 auto- 
correlation and cross-correlation galaxy and IA bias param- 
eters and 10 cosmological parameters, giving a total of 110 
parameters. 

The large number of parameters, especially when using 
the full bias model, favors the use of a Fisher matrix ap- 
proach. For N parameters, only N + 1 samples are required 
to estimate the parameter covariance matrix, dj = F i J , 
with, 

Fij=TjTj ^ Cov _ ld % { i) 



dpi 



dpj 



(48) 



where &(l) = {Cf MU , C\ ' 3 , C\ ' 3 , C/ j } are the set of 
observables across all multipole bins, and redshift bin com- 
binations. We consider correlations for 10 ^ £ ^ 3000 in 50 
logarithmically spaced bins in £ space. 

To calculate the partial derivatives, d$>(£) / dp, we take a 
2% reduction in each parameter with non-zero fiducial value, 
and an absolute step of —0.02 for all parameters whose fidu- 
cial value is zero. We checked that the results are insensitive 
to the exact size of the step size, obtaining consistent results 
with 1% and 3% step sizes. 

The covariance matrix Cov~^ between two observables, 
in multipole bin with mid- value £ and spanning £ m in(£) ^ 
£' < £max {£) , is given by 



C, 



+ c 



a 



WiXj 



C, 



w i x j 



f^fsky 
W i X j 



m ax(D 

=<»»(«) 



(49) 



We have modified the public ly available Cos moMC 
(|Lewis fc Bridldl2002r ) and CAMB ltewis et al.l (120001 ) codes 
to calculate the Fisher matrix and the correlation functions 
for the future survey specifications, in light of the IA and 
modified gravity models. 

Both relaxing the assumption that gravity is described 
by GR on cosmic scales and adding in systematic uncertain- 
ties increase the degrees of freedom fit by the prospective 



data and hence can degrade the quality of the cosmological 
information obtained. In figure [3] we show these dual effects 
on the 2D marginalised constraints for the dark energy pa- 
rameters. The constraints shown are for all data combined: 
CMB, galaxy position, lensing shear and cross-correlations, 
with the conservative bias marginalisation model using an 
A^fc = N z = 5 grid. Either including IAs, or allowing a mod- 
ification to gravity, separately has a roughly comparable ef- 
fect on reducing the constraining power on wo and w a , with 
IAs having a slightly larger impact. The inclusion of intrinsic 
alignments, while not significantly changing the degeneracy 
direction, does noticeably weaken the modified gravity pa- 
rameter constraints. 

We quantify the constraining power of the surveys us- 
ing the covariance matrix for the parameters , dj = F^ 1 . 
The diagonal elements of the covariance matrix give the 1- 
cr measurement uncertainty in each parameter, cr, = yCu- 
A 2 x 2 submatrix of a pair of parameters, C(pi,Pj), then 
gives the figure of merit (FoM) that includes the covariances 
between the parameters, 

-1/2 



FoM(pi, Pj ) = det[C(pi, Pj )Y 



(50) 



With this definition FoM = 1/<t 2 jj, where <r e // is the geo- 
metric mean of the principal axes of the 2-dimensional error 
ellipsoid. Note this differs by a factor of 1/6.177T, from an- 
other commonly quoted FoM, the area of the 95% confidence 
ellipsoid in the 2D marginalised space. 

We consider the Dark Energy Task Force figure of merit 
on the equation of state (EoS), FoM(wo,w a ) both in the 
absence of modifications to gravity, and also when the mod- 
ified gravity parameters and included and marginalised over. 
We quantify constraints on the modified gravity parameters 
by considering an equivalent 2x2 FoM, FoM (Qo, Qo(l + 
Ro)/2). 

An alternative way to measure dark energy constraints 
is to consider the combined constraints from equation of 
state and modified gravity together, through considering the 
FoM from a 4 x 4 covariance submatrix, 

FoM(comb) = det[C(w , w a , Q , Q (l + R )/2)}~ 1/4 . (51) 

IWang et all (|2010h applied an analogous statistic to com- 
pare constraints on a multi-parameter, model-independent 
dark energy figure of merit. This FoM accounts for all co- 
variances between the equation of state and modified grav- 
ity parameters, but at the same time entangles the EoS and 
MG constraints, with their different dependencies on mea- 
surements of the expansion history and growth of structure. 
With this definition, for n parameters, our FoM gives a 
measure of l/a e // 2 , the mean e rror in the n-dime nsional 
confidence ellipsoid. We note that I Wang et al.l |2010l ) define 
a slightly different FoM = l/a e // n , giving a figure of merit 
that scales with the volume of the n-dimensional space. 

For a 2 x 2 covariance matrix, M, for parameters x 
and y, one can calculate a correlation coefficient r corr — 
o~xy/o~xO~ y = y/1 — det(M) 1 0-2.0-2 . By dividing the 4x4 co- 
variance matrix for all 4 parameters, {wo, w a , Qo, Qo(l + 
Ro)/2}, into 2x2 submatrices, we can define an equivalent 
correlation coefficient between the equation of state param- 
eters and modified gravity parameters, 



1 - 



FoM(vjQ,w a )FoM(Qo, 



Qad+Rg) -, 
2 l 



FoM (comb) 2 



(52) 
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Figure 3. A comparison of dark energy constraints in the 2D marginalised parameter planes when intrinsic alignments are included 
using the LA model [full lines], in comparison to when it is assumed that IA's are perfectly understood and can be extracted to reveal 
the underlying cosmological shear and galaxy position correlations [dotted lines]. 95% confidence level constraints are shown when both 
GR is assumed [red] and when large scale modifications to gravity "MG" are allowed [black]. These results combine Planck- like CMB 
data with a Stage IV survey's galaxy position and shear auto and cross correlations. 



Table [3] summarises the 1-a, figure of merit (FoM) and 
correlation coefficient, r corr results for constraints coming 
from prospective Stage III and Stage IV surveys. We com- 
pare the different IA treatments, and the effect of lensing 
calibration and photometric redshift offsets. 

In the absence of modifications to gravity, IAs still have 
a significant impact, as was, for example, pointed out in 
(|joachimi fc Bridlell2010h . If IAs are assumed to be perfectly 
understood then one can achieve a ~ 14 fold improvement 
in the dark energy FoM from the Stage IV survey relative to 
Stage III. However when astrophysical uncertainties about 
IAs are included, and marginalised over, we find that the 
relative improvement of the photometric Stage IV survey, is 
reduced to 9. 

When the modification to gravity described in 12.11 is 
included, measurements of the growth of structure no longer 
purely constrain wo and w a . With IAs excluded, the EoS 
figure of merit with modified gravity allowed is weakened by 
50% relative to GR for a stage III survey, and by 70% for 
stage IV. 

The inclusion of IA uncertainties reduces both the EoS 



and MG figures of merit by roughly a factor of 2 relative to 
those when IAs are excluded from the analysis. Both with 
and without IAs included, the modified gravity FoM for a 
Stage IV survey is roughly a factor 6 improvement over that 
for Stage III. 

Overall, when both modified gravity and dark energy 
parameters are considered together, the FoM improves by a 
factor 8.5 between Stage III and Stage IV in the absence of 
IAs, and this is reduced to just under 7 with the conserva- 
tive modeling of IAs. The similarities in the adjustments in 
DE and MG FoM between Stage III and Stage IV suggest a 
high degree of correlation between the two parameter pairs. 
This is quantitatively reinforced by the correlation coeffi- 
cient r corr ; the correlations are higher for Stage IV than for 
Stage III and are degraded by ~15 and 25%, respectively, 
with the inclusion of intrinsic alignments. 

There is only a small difference between the figures of 
merit for the linear alignment and nonlinear alignment mod- 
els are included; with the NLA model giving slightly poorer 
constraints. The small difference suggests that the differ- 
ences between the LA and NLA models are to a large ex- 
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Survey 


Scenario 


<r(n A ) 


a(wo) 


cr{w a ) 


a(w p ) 


FoM(EoS) 


tr(Qo) 


^ /<9o(i+Ro)^ 


FoM(MG) 


FoM(comb) 


l^corr | 


Assuming GR 




no IA 


0.014 


0.169 


0.451 


0.029 


76.1 












Stage III 


LA 


0.025 


0.272 


0.677 


0.042 


35.2 














NLA 


0.027 


0.294 


0.728 


0.043 


32.3 














no IA 


0.003 


0.045 


0.127 


0.008 


1041.1 












Stage IV 


LA 


0.008 


0.081 


0.207 


0.017 


287.3 














NLA 


0.008 


0.086 


0.217 


0.017 


268.7 












Allowing an alternate modified gravity model 




no IA 


0.019 


0.226 


0.605 


0.034 


48.8 


0.666 


0.091 


40.1 


54.9 


0.59 


Stage III 


LA 


0.030 


0.363 


0.923 


0.043 


25.3 


0.894 


0.118 


25.0 


29.6 


0.53 




NLA 


0.032 


0.383 


0.971 


0.043 


23.8 


0.897 


0.118 


23.7 


27.7 


0.51 




no IA 


0.007 


0.064 


0.176 


0.017 


334.5 


0.198 


0.032 


253.6 


469.8 


0.79 


Stage IV 


LA 


0.010 


0.105 


0.274 


0.022 


166.1 


0.280 


0.041 


151.6 


202.1 


0.62 




NLA 


0.011 


0.106 


0.276 


0.022 


164.0 


0.292 


0.041 


143.7 


190.8 


0.60 




no IA 


0.008 


0.072 


0.204 


0.024 


200.5 


0.226 


0.039 


196.3 


327.9 


0.78 


Stage IV 


LA 


0.012 


0.111 


0.295 


0.031 


109.9 


0.290 


0.045 


129.8 


162.7 


0.64 


+sys. offsets 


NLA 


0.012 


0.112 


0.296 


0.031 


108.4 


0.301 


0.046 


123.8 


153.9 


0.61 



Table 3. Comparison of figures of merit and marginalised 1-cr errors for dark energy equation of state (EoS) parameters {wo,w a } and 
modified gravity (MG) parameters {QoiQo(l + Ro)/2}- A combined figure of merit including covariances between all 4 dark energy 
parameters, FoM (comb), and a correlation coefficient, r corr , between the EoS and MG parameters are also included. The table shows 
prospective constraints from galaxy position and weak lensing auto- and cross-correlations from Stage III and Stage IV surveys in 
combination with temperature and polarisation data from a Planck-like CMB survey. We assume a conservative model for galaxy and 
IA bias parameters, with A r j, = N z = 5. The results with "+sys. offsets" include marginalisation over weak lensing calibration and 
photometric rcdshift offset parameters meant to reflect possible instrumental systematic errors. 



tent drowned out by the uncertainties in the IA bias model. 
Small scale galaxy position correlations, which would be sen- 
sitive to differences in LA vs NLA through the gl term, are 
typically excluded since the multipoles exceed £ ma x- 

For the analysis shown in the table alone, we also con- 
sider the impact of two additional instrumental systematics: 
photometric redshift offsets and lensing shear calibration off- 
sets , on the figures of merit. Photometric redshift offsets, 
Azi alter the galaxy distribution inferred from observations 
as in (|16[) . Wh en systematic offsets a re considered, we model 
them following lAlbrecht et"al] (|2009l ): we allow independent 
offsets in each photometric redshift bins and impose a prior 
on these offsets of a(Azt) = 0.002. We model shear calibra- 
tion offsets by altering the measured shear correlations 

C e e '^'° ffset = (l + Am,)(l + Am,)C; ,e3 (53) 
C"^' o//se< = (1 + Am J )C; ,Ej (54) 

and impose a prior of a(Ann) = 0.001 N v h in each bin. 
Shear and redshift calibration offsets introduce an additional 
10 parameters. Shear calibration offsets are qualitatively de- 
generate with the inclusion of IA correlations of unknown 
amplitude in the ee correlation, and both cause a reduction 
in the figures of merit. While the systematic shear offsets 
in ee and ne are wholly correlated, the IA contributions can 
differ through the inclusion of the n degree of freedom in 
the GI cross-correlations. As such, the degradation in the 
constraints from including the instrumental systematic off- 
sets, as we model them here, are not as severe as those from 
marginalising over the uncertainties in the IA model. 

In figure [H we breakdown the impact of including IAs 
on the figures of merit derived as one combines CMB plus 
galaxy, lensing and galaxy lensing cross correlations in a 



piece- wise fashion. The inclusion of intrinsic alignments sig- 
nificantly deteriorates the expected dark energy constraints 
coming from weak lensing on its own, while the use of a grid 
bias model leads to the galaxy-galaxy correlations providing 
little constraining power on the non-bias parameters in the 
model. When lensing and galaxy position data are added in 
tandem, however, they are able to provide improved con- 
straints, over and above the lensing data alone. When IAs 
are included, the combined constraints are noticeably weaker 
than the constraints predicted by pure shear-shear measure- 
ments when IAs are neglected. The effects are mitigated to 
a good degree by gaining complementary information about 
the underlying cosmological potentials, and isolating out the 
IAs, by adding in galaxy position data. In particular, the in- 
clusion of cross-correlations, between galaxies and lensing, 
allow the correlated effects of the and GI IA contribu- 
tions to reduce the uncertainties in the IA model. 

Our findings for constraints on equation of state param- 
eters are consistent with those of IJoachimi fc Bridle! (|2010l) ; 
in the absence of CMB data, figures of merit with all galaxy 
position and weak lensing correlations and IAs included are 
comparable with those predicted by weak lensing alone in 
the absence of IAs. We do find, however that when we in- 
clude CMB data the FoM with IAs never becomes com- 
parable with those when IAs are excluded, even when all 
cross-correlation information is included. 

In a number of recent analyses of dark energy con- 
straints from prospective surveys, the uncertainties in the 
galaxy bias model are treated by a single, scale and red- 
shift independent, factor. In figure [5] we highlight that this 
assumption can have a dramatic effect on the predicted con- 
straining power of the survey. Allowing a gridded galaxy bias 
model, while excluding IA uncertainties, reduces the EoS fig- 
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Figure 4. How figures of merit (FoM) are affected by intrinsic alignments and the choice of data sets utilised in the analysis for 
equations of state parameters (wo vs w a ) and modified gravity parameters (Qo vs. Qo(l + Ro)/2)- We compare analyses in which IAs 
arc ignored [black full, triangle], where they are included using the linear alignment model [red dotted linos, square] and the non-linear 
alignment model [blue dashed, cross]. Datasets include a Planck-like CMB survey, denoted 'cmb', and Stage IV galaxy position-position 
'nn', shear-shear 'ee' and shear-galaxy position cross-correlations'ne'. 



ure of merit by roughly a third, and MG figure of merit by 
almost a half. 

Including intrinsic alignments in the analysis, while as- 
suming single, scale and redshift independent, amplitude has 
a marginal impact on the EoS and MG FoMs. In fact, inter- 
estingly, assuming that you know how IAs are formed and 
evolve provides additional information, actually improving 
the constraints. 

One can understand this by noting that the GG, II and 
GI components all depend on the underlying matter distri- 
bution but each exhibit a different evolution with redshift. 
Assuming a single-parameter normalisation, but multiply- 
binned measurements (from tomography) for the IAs, en- 
ables us to obtain an independent measurement of tp(zf) 
from the GG and both IAs. When uncertainties in the IA 
model are introduced, however, by using the grid bias model, 
they markedly degrade the dark energy constraints. One 
doesn't see the plateauing of the FoM that one would see 
with intrinsic alignments excluded from the analysis. 

We rationalise the plateau with no intrinsic alignments 
as follows: the contributions from shear-shear, shear-position 



and position-position power spectra have different depen- 
dencies on redshift resolution. We expect the constraining 
power of position-position alone to be weak given the multi- 
pole cuts, and to be badly affected by even a small number of 
free bias parameters. Therefore we should be dominated by 
shear-shear and shear-position information which are much 
more resilient to a lack of redshift information because, when 
IAs are not included, each is modulated by the broad lensing 
weight function. 

We can combine the information in figures 0] and [S] to 
identify how many free parameters can be accommodated 
before constraints degrade relative to conventional con- 
straints: from shear-shear correlations alone, in which IAs 
are ignored. CMB plus shear-shear data alone give a modi- 
fied gravity figure of merit of 170 (figure[4]), which is roughly 
the same value obtained when including all two-point cross- 
correlations and including a bias grid with iVbias = 4 (figure 
[5|. For the equation of state figure of merit we can use a 
bias grid with Afbias = 3 before we reach the same figure of 
merit using all two-point functions as we would obtain from 
the traditional approach. A bias grid with Nbi& s = 3 has 
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Figure 5. The impact of the number of k and z bins, = iVz = N^i aa , in the bias model on the equation of state (EoS) and modified 
gravity (MG) figures of merit (FoM). Scenarios are shown in which IAs arc excluded [black .triangle], and in which linear alignment 
(LA) [red, square] and nonlinear alignment (NLA) [blue, cross] models for intrinsic alignments are used. If IAs are excluded one sees 
a platcauing of the figure of merit as the number of bias marginalisation parameters in increased. With the addition of parameters 
to describe uncertainties in the IA amplitude no such plateauing is seen. The inclusion of IA, with an assumption that they are well 
understood, and can be described by scale and redshift independent nuisance parameters (JV(,j os = 1) actually improves the dark energy 
constraints because the IAs provide additional cosmological information about the high redshift potential <f>(zf). If uncertainties in the 
IA model are allowed however, there is a significant deterioration in the constraints on both FoM. The results presented here are for 
prospective CMB and Stage IV large scale structure survey utilising all galaxy position and shear auto- and cross- correlations. 



3 x 3 x 4 = 36 free parameters in total, 18 for galaxy bias 
and a further 18 for the intrinsic alignment model. We there- 
fore need astrophysics to be sufficiently kind that the bias 
functions are sufficiently smooth in both scale and redshift, 
or to have sufficient information from simulations to be able 
to parameterise the functions with roughly this number of 
free parameters. 

The inclusion of information in the mildly nonlinear 
regime can have a potentially significant effect on improving 
constraints on dark energy parameters, purely as a result of 
the large number of modes available to include in the anal- 
ysis. If modifications to gravity are included however then 
one has to make an assessment of how well large scale struc- 
ture growth in the non-linear regime is understood. Recent 
analyses show that in some specific modified gravity theo- 
ries there can be subtleties in the non-linear behaviour that 



might have to be included JOyaizu 2008 


; Ovaizu et al. 20081; 


Schmidt et al.l 20091; 


Khourv & Wvman 


20091; Ferraro et al.1 


20111; ICui et al.ll201C 


;IZhao et al.ll201ll;lBrax et alj|201ll). In 



figure [6] we highlight the sensitivity of the figure of merit to 
the assumptions about the smallest scales to be included in 
the analysis, parameterised here by l ma x- 



When CMB data is included, the pressure to go to high 
multipoles is reduced. There is only a 50% increase in FoM 
on increasing the maximum multipole from 1000 to 3000, 
compared to over a factor of two when CMB data is not 
used. 



4 CONCLUSIONS 

Using tests of the expansion history of spacetime and the 
growth of large scale structure, in tandem, gives the best 
prospects for testing gravity on cosmic scales. Weak lensing, 



galaxy position, CMB ISW and potential peculiar velocity 
observations provide very complementary constraints on the 
gravitational potentials, through measuring both their sum, 
4> + tp and ip on its own. 

Fundamental to realising the full potential of these 
complementary observations is a requirement to minimise 
both instrumental and astrophysical systematic uncertain- 
ties that can dilute the cosmological constraining power of 
upcoming surveys. Weak lensing observations could poten- 
tially offer a direct way to measure the gravitational poten- 
tials without the bias uncertainty in relating galaxy posi- 
tions to the underlying CDM matter distribution. On the 
other hand, intrinsic alignments provide a significant sys- 
tematic signal. Uncertainties about IA formation, evolution, 
and variation amongst galaxy-type, have to be factored into 
a realistic assessment of how well weak lensing shear mea- 
surements can constrain a cosmological model. 

In this paper, we have shown that how systematic un- 
certainties are modelled can have a profound impact on the 
predicted dark energy and modified gravity constraints from 
future large scale structure imaging surveys. 

By utilising a grid-based approach to marginalise over 
uncertainties in both galaxy bias and intrinsic alignment 
contributions to lensing shear and galaxy position correla- 
tions, we have provided conservative and optimistic bounds 
for constraints on both dark energy equation of state param- 
eters and a useful phenomenological modified gravity model. 

We considered three figures of merit to quantitatively 
compare constraints on the dark energy equation of state 
and modified gravity parameters both separately and in 
combination. Quoting separate figures of merit for EoS and 
MG parameters can be used to show their different depen- 
dences on data sets and assumptions. They can also contrast 
the equation of state parameter dependence on expansion 
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Figure 6. The impact of including observations on small scales, denoted by the maximum multipole, l m ax, up to which correlations 
are considered, on the equation of state (EoS) and modified gravity (MG) figures of merit (FoM) . Results are shown for a Stage IV 
photometric survey alone [upper panel] and [lower panel] including complementary constraints from a Planck-like CMB survey when 
IAs are excluded [black full line] and included using the LA [red, dotted] and NLA [blue, dashed] models. While including smaller-scale 
observations would appear to improve both figures of merit, one has to consider the theoretical uncertainties present in modeling these 
small scales, especially in the context of modifications to gravity, therefore it is worthwhile assessing how a conservative approach of 
neglecting such scales might impact the projected cosmological constraints. 



history measurements, and modified gravity parameter de- 
pendence on the growth history. We have found, however, 
that there is a high degree of correlation between these two 
sets of parameters so that treating them independently ig- 
nores an important association. We proposed, and quanti- 
fied, a combined figure of merit and related correlation co- 
efficient as a way to address this. 

We have found that the constraints have a significant 
sensitivity to how galaxy bias and intrinsic alignments are 
incorporated into the analysis. The equation of state and 
modified gravity figures of merit are a factor of 4 smaller 
when a conservative scale and redshift dependent grid model 
is used, than when bias and IAs uncertainities are assumed 
to be redshift and scale independent. Marginalising over sys- 
tematic uncertainties in the IA model led to a factor of two 
reduction in the figures of merit. 

Whether a linear alignment or nonlinear alignment 
model underpinned the IA model had only a small ef- 
fect in comparison to our assumptions about the evolu- 



tion of bias and IAs in redshift and scale. Understand- 
ing the astrophysical evolution and population dependence 
of intrinsic alignments, therefore, could dramatically im- 
prove the cosmological information that comes out of fu- 
ture photometric large scale structure surveys, such as DES, 
SuMIRe, Euclid, LSST and WFIRST. We discuss the impli- 
cations of weak lensing systematics for optimising cosmic 
shear surveys to measure dark energy in a rel ated paper 
(j D. Kirk. I. Laszlo. S. Bridle and R. Beanll20lil ). 

In addition to uncertainties in bias and IAs, an under- 
standing of evolution in the nonlinear regime could also have 
a profound impact on constraints, through increasing the 
maximum multipole to which analyses can proceed. If GR 
governs cosmic evolution this may be achievable, while the 
model dependence of the non-linear regime in modified grav- 
ity models could make this far more challenging. 

Combining information from the photometric surveys 
we have considered here with that from spectroscopic galaxy 
data, such as might come from BigBOSS, EUCLID and 
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WFIRST, might allow closely situated galaxies to be iso- 
lated, and their intrinsic alignments to be studied. We will 
consider in future work how this could provide an important 
avenue to improve our understanding of intrinsic alignments 
and in turn maximise the cosmological constraining power 
of future wide and deep photometric surveys. 

We finally note that, as part of this work, we have pro- 
vided a fitting function to allow other researchers to generate 
weak lensing and galaxy position correlations for modified 
gravity theories of the form we consider here. 
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A APPENDIX: MODIFIED GRAVITY 
FITTING FUNCTION 

A.l Fitting function form 

Some modified gravity (MG) models, such as f(R) theo- 
ries, can be tailored to reproduce a selected expansion his- 
tory. However their predictions for the growth of structure 
can then differ from that predicted by that expansion his- 
tory assuming GR. Here we consider a fitting function for 
a modified gravity model in which the expansion history 
is described by ACDM, but the growth history is modified 
through a deviation from GR described by two parameters, 
Qo and Ro, and a third parameter s which encapsulates the 
time dependence of the deviation: 



ACDM for the same cosmological parameters: 



Q(a) = f + (Q -l)a s 
R(a) = l + {Ro- l)a s 



(55) 



In the main analysis in this paper we have assumed s = 3. 

The key input into calculating the galaxy position and 
weak lensing shear observables is the matter power spectrum 
as a function of scale, k, and redshift, z. Here we obtain a 
analytical fit for the ratio of the matter power spectrum 
in the modified gravity model in (|55[) . to that predicted by 



rfit(Qo, Ro,s) = 



P(k) t 



(56) 



in, AC D M 

To motivate the form of the fit, we note that 
the behaviour of growth in this model is described in 
I Bean fc Tangmatithaml < |201CT > has two distinct regimes as 
given in ([7]) and JSJ. On small scales the growth purely de- 
pends on the product QoRo, via 0. On large scales, the be- 
haviour involves various derivatives of the modified gravity 
parameters, and the equations depend uniquely on Qo and 
Ro- Our fit distinguishes between these two regimes in scale, 
and fits the evolution with redshift of the high wavenumber 
C-fT) and low wavenumber ('I/') regime separately. The two 
regimes are joined via a third function, x(k), assuming a 
transition scale, k c - 



rfu(k, z) 
r L (z) 



r H {z) 
x(k) 



[l-ir(fc)] r L (z) + x(k)r H (z) 
1+ \L 1 (1-R )+L 2 (e 1 



')] 



1 + H^l-QoRo) [e" 2Z + H 3 
k \ v 

tanh 



(57) 

(58) 
(59) 
(60) 

(61) 



The values for the 10 fitting parameters 
{Li, 1/2, L-j,, 1/4, Lb, Hi, H2, H3, k c ,p}, are obtained us- 
ing OriginLab's Origin software to fit our custom function 
to the power spectra coming from the CAMB code. While 
s = 1 and s = 3 have been most commonly used choice in 
the literature, we obtain the fit for s — 1 — 4, using a grid of 
values in Qo and Ro between 0.9 and 1.1. The spectra were 
calculated at 50 redshift steps in ^ z ^ 3 and for over a 
hundred values in k ranging from 7 x 10~ 6 to 40Mpc _1 . 
For s = 3, 23 different {Qo,^o} pairs were used to obtain 
the fit, for s = 1 and s = 2 a subset of these was used, and 
found to be sufficient to achieve sub percent accuracy in 
the Cis. 

Table |A~T1 provides the values for the fits parameters for 
each value of s. 



A. 2 Performance 

The fit given above reproduces the ratios of the matter power 
spectrum in (|56|) . derived from CAMB, to better than 0.01% 
for k > lOfcc, within 0.6% around the transition scale, k c , 
and within 0.8% at scales around the horizon scale today. 
This is sufficient to reproduce the C(S in the modified gravity 
scenario to within sub percent (~0.1%) levels. The error is 
largest for low la and increases as s decreases. 

Figure [7] shows the ratio of the matter power spectrum 
today for the modified gravity model to that for the fiducial 
ACDM model for both the fit and the full integration using 
CAMB. At z = 0, where the modifications are largest, as a 
function of k, the fit matches the simulations to within 0.8 
% and remains accurate at this level for all 2 for 5% changes 
in Qo and Ro- 

When calculating the angular correlation function, Ci, 
the power spectrum fit is integrated over k, z. Applying the 
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Fit Parameter 


s = 1 


s = 2 


s = 3 


s = 4 


Lil 


0.5293 


n a n a T 
0.494/ 


0.4268 


0.3635 


£-2 


-4.733 


-5.692 


-6.300 


-6.575 


U\ 


-1.610 


-1.660 


-1.764 


-1.817 


L4 


-0.8678 


-1.610 


-2.409 


-3.263 


L5 


0.2878 


0.0867 


0.0336 


0.0156 


Hi 


-0.5655 


-0.2023 


-0.0984 


-0.0557 


H 2 


-0.6144 


-1.263 


-1.935 


-2.718 


H 3 


0.1803 


0.0754 


0.0317 


0.0153 


kdMpc- 1 ] 


1 X 10~ 3 


7 X 10-' 


1 6 x 10~ 4 


5 x 10~ 4 


P 


0.9422 


1.048 


1.076 


1.090 


Table 4. Summary of values for the 10 parameters used in the 
fitting function, given for each value of s, the power law exponent 
in the modified gravity function 155t. 
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Figure 7. Ratios of z = matter power spectra in modified 
gravity model to the fiducial ACDM model obtained via simu- 
lation [red] compared with the ratios obtained using the fitting 
function [blue]. Two models shown are Qo = 1, Ro = 0.95 [dashed 
lines] and Qo = 0.95, Ro = 1 [dotted lines]. At small scales the 
two models are degenerate, since the evolution of the matter per- 
turbations is only dependent on the product QR, while at large 
scales their behaviours arc distinct. The fitting function provides 
agreement to within 0.01% for most scales. At the transition scale 
k ~ k c and on horizon scales the fit is a litter poorer, ~ 0.8%, how- 
ever this limited range of scales contributes only a small amount 
to the angular correlations used in the analysis. 



fit factor to a standard power spectrum, and using this 
fit power spectrum to compute Cts for galaxy autocorrela- 
tions, galaxy-weak lensing cross-correlations, and weak lens- 
ing auto-correlations results in sub percent(0.1%) level ac- 
curacy, as shown in figure [8] 
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